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We consider the classical magnetoresistance of a Weyl metal in which the electron Fermi surface 
possess nonzero fluxes of the Berry curvature. Such a system may exhibit large negative magneto- 
resistance with unusual anisotropy as a function of the angle between the electric and magnetic 
fields. In this case the system can support a new type of plasma waves. These phenomena are 
consequences of chiral anomaly in electron transport theory. 
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Materials with nontrivial topological properties have 
attracted considerable interest after the discovery of 
topological insulators [1]. One type of such materials is 
the so-called Weyl semimetals, characterized by the the 
presence of points of band touching (Dirac points) [2- 
10]. In this paper, we study the metallic counterparts of 
these materials — the Weyl metals, where Dirac points are 
hidden inside a Fermi surface. We show that these mate- 
rials may exhibit large negative magnetoresistance with 
unusual anisotropy. We also find a new type of plasma 
waves in these systems. 

At low magnetic field B and at relatively high temper- 
ature Landau quantization can be neglected, and electron 
transport in metals can be described using the semiclas- 
sical Boltzmann kinetic equation 
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Here n p (r, t) is the electron distribution function, p is 
the quasimomentum, I C oii{n p } is the collision integral, 
and 



f=^+pxn p , 



eE 



-f x B. 



(2a) 
(2b) 



The last, "anomalous," term in Eq. (2a), proportional to 
the Berry curvature 
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was introduced in Rcf. [11]. (See also reviews on the 
subject in Refs. [12, 13].) In systems with time-reversal 
symmetry f2 p = fi_ p , while in centro-symmetric sys- 



tems f2 t 



-f2_ p . Thus, in systems which are both 



time- and centro-symmetric fl p = 0. In this case the 
magneto-resistance described by Eq. (1) is positive and 
is governed by the parameter (u; c Tt r ) 2 [14]. Here cu c is 
the cyclotron frequency and r tr is the electron transport 
mean free path. The Berry curvature is divergence-free 
except at isolated points in p space, which are associated 
with band degeneracies. As a result, in the case where 



the electronic spectrum has several valleys, they can be 
characterized by integers (see, for example, Ref. [15]) 
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Here the index "i" labels the valleys, dS is the elementary 
area vector. Nonzero values of fcW are realized if near 
the degeneracy points, electrons can be described by the 
massless Dirac Hamiltonian [27] 



H = ±V(T ■ P. 
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Here P = — fiV — |A is the momentum operator, A is 
the vector potential, er is the operator of pseudospin, v 
is the quasiparticle velocity, and the signs ± correspond 
to the different chiralities of the Weyl fermions. 

It is well known that massless Dirac fermions ex- 
hibit chiral anomaly which can be understood in the 
language of level crossing in the presence of a mag- 
netic field [16]. According to the Nielsen-Ninomiya theo- 
rem [17], the number of valleys with opposite chiralities 
(positive and negative values of fcW) should be equal, 
and so ^ feW = 0. Recently, gapless semiconductors 
with topologically protected Dirac points (Weyl semi- 
metals) have attracted significant attention [2-10]. Both 
time-reversal-breaking [2], and non-centro-symmctric [4] 
versions of these systems have been proposed. In the 
absence of a random potential and doping, the chemical 
potentials in these systems are at the Dirac points. 

In this article we consider the case where in equilib- 
rium the chemical potential fi = fa measured from the 
Dirac points is finite, and show that the semiclassical 
Eqs. (1) and (2) can yield a substantial anomaly-related 
negative magnetoresistance. The latter also exhibits un- 
usual anisotropy as a function of angle 9 between E, and 
B. Here fa is the chemical potential in the i-th valley, 
measured from the Weyl's point. Using Eqs. (2) we get 
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where v = de p /dp. Substituting Eqs. (6) into Eq. (1) we get the kinetic equation in the form 
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(cf. Ref. [18]). Let us consider the case where /i » T, 
ftw c — h\e\v 2 B/cn, and assume that the conductivity of 
the system is determined by elastic scattering. Then, the 
collision integral I co ii in Eq. (7) describes elastic intra- 
and inter-valley scattering. We assume that r tr <C r, 
where r is an elastic inter-valley scattering mean free 
time. In this case the anisotropy of the intra- valley distri- 
bution function can be neglected, and the latter depends 
only on the energy e: n$ — nW(e). Denoting by p^(e) 
the density of states [12, 13], 

in the homogeneous case we get the kinetic equation in 
a form 

dn^(e) jfeW e 
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where the collision integral now includes only inter-valley 
scattering. For this, we will use the relaxation time ap- 
proximation 
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where dn l (e) is the deviation of the distribution function 
from its equilibrium value. 

The electron density and entropy density in i-th valley 



are 
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Integrating Eq. (9) over pW(e)de we get the conservation 
law for particle number in each valley: 

(E-B)-—- , (13) 
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Thus, in the presence of electric and magnetic field, the 
number of particles in the i-th valley, JVW, is not con- 
served even if r — > oo. This is the chiral anomaly which 
was originally introduced in field theory in Refs. [19], and 
latter discussed in the context of electron band struc- 
ture theory in Ref. [16], and in the theory of supcrfluid 
3 He [20, 21]. It is interesting that the anomaly can be 
understood completely in the framework of semiclassical 
kinetic equation Eq. (1), characterized by fcW [22], and 
that the term proportional to E • B in Eq. (13) is the 
same as obtained in Ref. [16] in the ultra-quantum limit. 

Equation (7) represents a low energy effective theory. 
To see why the number of electrons in an individual val- 
ley is not conserved one has to take into account the 
spectral flow process which bring the energy levels (to- 
gether with electrons occupying them) from one Dirac 
point to another through the bulk of the valence band, 
as schematically shown in Fig. 1. Such a possibility exists 
only in the presence of a magnetic field. 

The existence of the chiral anomaly results in a rather 
unusual mechanism for negative magnetoresistance. The 
easiest way to calculate the magnitude of the effect is to 
estimate the rate of entropy production in the presence 
of electric field. 
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At small E, the stationary solution to Eq. (9) is 
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For simplicity let us assume there are only two valleys, 
with fci ; 2 = ±1, and with the same quasiparticle velocity 
v. Let the z-axis be parallel to B. Then, from Eqs. (15) 
and (16) we get an anomaly-related contribution to the 
component a zz of the conductivity tensor 
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Note that a zz given by Eq. (17) is an increasing function 
of the magnetic field. All other anomaly- related compo- 
nents of the conductivity tensor cr^ are zero. In other 
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words, the anomaly-related current can flow only in the 
direction of B. One can also understand this fact by 
noticing that, at E = 0, and in the presence of a mag- 
netic field, Eq. (14) gives an expression for a current den- 
sity [23, 24] 
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Here we assume electron distribution functions in the in- 
dividual valleys have equilibrium forms. In the case of a 
global equilibrium, all /Ltj = /i, and the contributions to 
Eq. (18) from different valleys cancel each other. Accord- 
ing to Eq. (13), in the presence of electric and magnetic 
fields, an imbalance of electron populations and, conse- 
quently, a difference between the the chemical potentials 
[Mi is created. As a result, there is a finite current den- 
sity, which can relax only via intra-valley scattering. In 
agreement with Eq. (17), its value is proportional to r, 
its direction is parallel to B, and it responds only on the 
component of the electric field parallel to B. 

There is a significant difference between the anomaly- 
related [Eq. (17)] and the conventional Drudc contribu- 
tions cr^(B) to the B-dependence of the conductivity 
tensor. For an isotropic Fermi surface and in the relax- 
ation time approximation, all components of <jf^\ except 

for <riz\ are decreasing functions of B. For anisotropic 
Fermi surface, at (w c r tr ) 2 <C 1, there is a B-dependence 



of a zz as well, which can be estimated as 
^)(0)- ff (?(B)~^(0)( Wc r tr 
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Herecr( D )(0) = e 2 ^v 2 Tt r /3 is the Drude conductivity, and 
v ~ /v 3 is the density of states at the Fermi level. 
Thus, for small magnetic field, the anomaly-related con- 
tribution (Eq. 17) dominates the magnetoresistance, pro- 
vided 
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Generically, in small gap semiconductors, the parameter 
T~/T tr 3> 1, because the inter- valley scattering requires 
a large momentum transfer. If the scattering potential 
is smooth, this parameter become exponentially large. 
Even in the case of an anisotropic Fermi surface, de- 
pending on symmetry there could be a direction of E 
for which a'fz is independent of B. At small values of 
yu <C T the conductivity is determined by electron-hole 
scattering [6, 7]. In this case one has to substitute ji for 
T in Eqs. (17) and (20), while the parameter r/T tr ^> 1 
is exponentially big. 

For ((x> c T tr ) 2 3> 1 the B-dependence of cr^(B) sat- 
urates and it becomes independent of B. In contrast, 
the deviation of the anomaly-related contribution from 
the quadratic in B behavior takes place at much higher 
magnetic fields. Thus, a zz could be a non-monotonic 



function of B. Finally, the anomaly-related contribution 
to the conductivity tensor may be distinguished by its 
unusual frequency dependence: it is controlled by the 
parameter (lot) 2 , rather than by the conventional pa- 
rameter (uiT tr ) 2 . Here u> is the frequency of the electric 
field. 

At low values of fx the anomaly-related contribution 
to the conductivity can be even bigger than the Drude 
contribution a^. In this case system supports a new 
type of weakly damped plasma waves with a frequency 
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provided cu p 3> t _1 . 

The approach based on the semiclassical equations of 
motion, Eqs. (1) and (2), is valid if ll 3> hw c . In the 
opposite, ultra-quantum, limit w c Tt r ^S> 1, the anomaly- 
related negative magnetoresistance has been previously 
discussed in Refs. [16, 25]. In this case the spectrum of 
the Dirac equation has the form 
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For n — case eo = ± w Pz, where ± corresponds to dif- 
ferent valleys. In other words, the n — Landau level 
is chiral: the branches of the spectrum with e = ±vp z 
correspond to different valleys, as shown in Fig. 2. Con- 
sider the case where both the chemical potential and 
the temperature are small compared to the energy dif- 
ference between the zero and the first Landau levels, i.e. 
H, T < Hv/Lb, where Lb = y/hc/eB is the magnetic 
length. In this case only chiral branches of the spectrum 
are occupied by electrons. Contributions to the current 
from branches of the spectrum with different chiralities 
can relax only by inter- valley scattering processes charac- 
terized by r. If the electric field is applied in the z direc- 
tion, electrons movie according the law p z = eE z —p z /r; 
v z = ±v, and we get the following expression for the 
conductivity [16]: 



(23) 
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By the same token we can obtain an expression for the 
plasma frequency at zero wave vector, 
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where r s = e 2 /nhv, and k is the dielectric constant. 
Eq. (24) is valid if r s < 1. 

The fact that Eq. (24) does not have a classical limit 
(h — > 0) is a particular example of a general property of 
collective modes in the massless Dirac plasma [26]. We 
note however, that in three dimensions and at B = 
the plasma frequency is proportional to fi. In contrast, 
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FIG. 1: Schematic 3D electron spectrum in a Weyl metal. 
Only two valleys in the electron spectrum are shown. The 
dashed line indicates the direction of the electron spectral 
flow in the presence of parallel electric and magnetic fields. 




FIG. 2: Schematic electron spectrum of a Weyl metal in the 
ultra-quantum limit, n = 0, 1 label Landau levels. The 
dashed line indicates the direction of the electron spectral 
flow in p z space in the presence a z-component of the electric 
field. 



Eq. (24) is independent of the value of /i, and remains 
finite even when fj, = 0. 

Equations (23) and (24) are valid only for Weyl metals 
where time reversal symmetry is preserved. In systems 
with no time reversal symmetry linear in B contributions 
to <t zz are allowed. The magnitude of these contributions 
is not universal and depends on the details of the mech- 



anism of time-reversal symmetry violation. 
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